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Study of the decay B°(B°) — > p + p , and constraints on the CKM angle a. 
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Using a data sample of 89 million T(4S I ) — ► BB decays collected with the BABAR detector at the 
PEP-II asymmetric B Factory at SLAC, we measure the B°(B°) — > p + p~ branching fraction as (30± 
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4(stat)±5(syst)) x 10 6 and a longitudinal polarization fraction of fh = 0.99±0.03(stat) to'o3( svs t)- 
We measure the time-dependent-asymmetry parameters of the longitudinally polarized component 
of this decay as C L = -0.17 ± 0.27(stat) ± 0.14(syst) and S L = -0.42 ± 0.42(stat) ± 0.14(syst). We 
exclude values of a between 19° and 71° (90% C.L.). 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The recently observed 0] decay B Q (B°) -> p+p- pro- 
ceeds mainly through the b — > uud tree diagram. Interfer- 
ence between direct decay and decay after B°-B° mixing 
results in a time-dependent decay-rate asymmetry be- 
tween B a and B° that is sensitive to the CKM angle 
a = axg[-V td V t t/V ud V* b ]. The presence of loop (pen- 
guin) contributions introduces additional phases that can 
shift the experimentally measurable parameter a e e away 
from the value of a. In the presence of penguin contri- 
butions a c ff — a + <5a pe nguin- A constraint on a tests 
the Standard Model description of CP violation. Recent 
measurements of the B + — > p~*_P° branching fraction and 
upper limit for B° — > 3] indicate small penguin 

contributions in B — > pp, as has been found in some cal- 
culations . Here we present a time-dependent analysis 
of B°(B a ) — > p + p~ . 

The CP analysis of B decays to p + p~ is complicated by 
the presence of three helicity states (h = 0, ±1). The h = 
state corresponds to longitudinal polarization and is 
CP-even, while neither the h = +1 nor the h = — 1 state 
is an eigenstate of CP. The longitudinal polarization 
fraction Jl is defined as the fraction of the helicity zero 
state in the decay. The angular distribution is 

d 2 T _ 
Td cos Videos #2 

\ [f L cos 2 9 1 cos 2 2 + 1(1 - f L ) sin 2 1 sin 2 2 )(1) 

where 6i,i — 1 , 2 is defined for each p meson as the angle 
between the ir° momentum in the p rest frame and the 
flight direction of the B° in this frame. We have inte- 
grated over the angle between the p-decay planes. A full 
angular analysis of the decays is needed in order to sep- 
arate the definite CP contributions; if however a single 
CP channel dominates the decay, this is not necessary [5j . 
The longitudinal polarization dominates this decay [ilia]- 
This measurement is based on 89 million Y(4S) — * BB 
decays collected with the BABAR Q detector at the 
PEP-II asymmetric B Factory at SLAC. We reconstruct 
B°(B°) — > p + p~ candidates (B rec ) from combinations of 
two charged tracks and two n° candidates. We require 
that both tracks have particle identification information 
inconsistent with the electron, kaon, and proton hypothe- 
ses. The 7T° candidates are formed from pairs of pho- 
tons that have measured energies greater than 50 MeV. 
The reconstructed n° mass must satisfy 0.10 < to 77 < 
0.16 GeV/c 2 . The mass of the p candidates, m^i^o, 
must satisfy \m n ±„o - 0.770 GeV/c 2 | < 0.375 GeV/c 2 . 
When multiple B candidates can be formed we select 



the one that minimizes the sum of the deviations of 
the reconstructed 7r° masses from the true 7r° mass. 
Combinatorial backgrounds dominate near | cos 6i\ = 1 , 
while backgrounds from B decays, like B° — > p + 7r~, 
with an additional low energy ir from the rest of the 
event (ROE), tend to concentrate at negative values of 
cos#i. We reduce these backgrounds with the require- 
ment -0.8 < cos6»i < 0.98. 

Continuum e + e~ — > qq (q = u,d,s,c) events are the 
dominant background. To discriminate signal from con- 
tinuum we use a neural network (Af) to combine six vari- 
ables: the two event-shape variables used in the Fisher 
discriminant of Ref . || ; the cosine of the angle between 
the direction of the B and the collision axis (z) in the 
center-of-mass (CM) frame; the cosine of the angle be- 
tween the B thrust axis and the z axis; the cosine of the 
angle between the B thrust axis and the thrust axis of 
the ROE; the decay angle of the 7r° (defined in analogy 
to the p decay angle, the sum of transverse momenta 
in the ROE relative to the z axis. 

Signal events are identified kinematically using two 
variables, the difference AE between the CM energy of 
the B candidate and y/s/2, and the beam-energy substi- 
tuted mass mEs = yj (s/2 + p; • ps) 2 / Ef — p 2 B , where 
y/s is the total CM energy. The B momentum ps and 
four-momentum of the initial state (-Ej,Pi) are defined 
in the laboratory frame. We accept candidates that sat- 
isfy 5.21 < m ES < 5.29 GeV/c 2 and -0.12 < AE < 
0.15 GeV. The asymmetric AE window suppresses back- 
ground from higher-multiplicity B decays. 

To study the time-dependent asymmetry one needs to 
measure the proper time difference, At, between the two 
B decays in the event, and to determine the flavor tag of 
the other i?-meson. The time difference between the de- 
cays of the two neutral B mesons in the event (B rec , Stag) 
is calculated from the measured separation Az between 
the B ICC and B tag decay vertices 0, U3 • We determine 
the -B roc vertex from the two charged-pion tracks in its 
decay. The Stag decay vertex is obtained by fitting the 
other tracks in the event, with constraints from the B rec 
momentum and the beam-spot location. The RMS reso- 
lution on At is 1.1 ps. We only use events for which the 
proper time difference between the B rcc and -Btag decays 
satisfies |At| < 20 ps and the error on At, er(At), is less 
than 2.5 ps. The flavor of the Bt ag meson is determined 
with a multivariate technique |8| that has a total effec- 
tive tagging efficiency of (28.4 ± 0.7)%. The events are 
assigned to five mutually exclusive taggging categories 
Lepton, Kaon 1, Kaon 2, Inclusive, and Untagged, 
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listed in order of decreasing reliability of the tag. 

Signal candidates may pass the selection even if one 
or more of the pions assigned to the p + p~ state is from 
the other B in the event. These self-cross-feed (SCF) 
candidates comprise 39% (16%) of the accepted signal 
for fh = 1 (h = 0). 

The efficiency of the selection is 7.7% (14.9%) for lon- 
gitudinally (transversely) polarized signal as determined 
with Monte Carlo (MC) The signal efficiency tak- 

ing into account the measured polarisation is 7.7%. We 
select 24288 events, which are dominated by combina- 
toric backgrounds: roughly 86% from qq and 13% from 
BB. We distinguish the following candidate types: (i) 
correctly reconstructed signal, for B°(B°) — > p + p~ de- 
cays where the correct particles are combined to form 
the B ICC candidate; (ii) SCF signal; (iii) charm B ± back- 
ground (b — * c); (iv) charm B° background (b — ► c); (v) 
charmless B backgrounds; (vi) continuum e + e~ — * qq 
(q = u,d,s,c) background. We consider both types (i) 
and (ii) as signal. The charmless decays J5 ± — * p^n , 
B ± — > p ± p°, B ± —>■ aiiT , and B ± '° — > aip are expected 
to contribute to the final sample. For these decays we 
assume the following branching fractions: 11.0 ± 2.7 [l2| . 
26A±l s 4 [13, 30 ±15, and 20 ±20, in units of 10" 6 , corre- 
sponding to 17 ± 4, 16 ± 4, 30 ± 15, and 26 ± 26 events in 
the data, respectively. The latter two are estimated from 
the measured branching fractions of related decays. We 
expect an additional 283 ± 283 candidates of charmless 
B decays with more than four mesons in the final state; 
since most branching fractions for such modes have not 
been measured yet, we generate them using the JETSET 
simulation [J^. We expect 1700 (1016) charged (neutral) 
B decays to final states containing charm mesons. The 
rest of the background is composed of continuum. Each 
of these backgrounds is included as an individual compo- 
nent in the fit, where the continuum yield is allowed to 
vary in the fit. 

Each candidate is described with the eight B rec kine- 
matic variables toes and AE, the m^i^o and cos 9i val- 
ues of the two p mesons, At, and AT. For each differ- 
ent candidate-type considered, we construct a probabil- 
ity density function (PDF) that is the product of PDFs 
in each of these variables, assuming that they are uncor- 
rected. The total PDF is used in the fit to data. 

The parameters of the PDFs for continuum- 
background toes j AE, cos 9i, and Afaie allowed to vary 
in the final fit to the data. The distribution of the con- 
tinuum as a function of tyi- k ± 7T o is described by a non- 
parametric PDF ^3 derived from toes and AE data 
sidebands. For all other types these distributions are 
extracted from high-statistics MC samples. The cos^ 
distributions for the background are described by a non- 
parametric PDF derived from the MC, as the detec- 
tor acceptance and selection criteria modify the known 
vector-meson decay distribution. The signal distribution 
is given by Eq. (1) multiplied by an acceptance function 



determined from signal MC. We take into account known 
differences between data and the MC. 

The signal At distribution is described by an exponen- 
tial (B lifetime) multiplied by a CP violating term, con- 
voluted with three Gaussians (~ 90% core, ~ 9% tail, 
~ 1% outliers) and takes into account a (At) from the 
vertex fit. The resolution is parameterized usinga large 
sample of fully reconstructed hadronic B decays |9j . The 
nominal At distribution for the B backgrounds is a non- 
parametric representation of the MC; in the study of sys- 
tematic errors we replace this model with the one used 
for signal. The continuum background is described by 
the sum of three Gaussian distributions whose parame- 
ters are determined by fitting the data. 

The signal decay-rate distribution /+(/-) for B tag = 
B° (B°) is given by: 

e -|At|/r 

f±(At) = [1 ± Ssm(Am d At) =f C cos( Ato^ At)] , 

4r 

where r is the mean B° lifetime, Am d is the B°-B° mix- 
ing frequency, and S— Sl or St and C= Cl or Ct are 
the CP asymmetry parameters for the longitudinal and 
transversely polarized signal. The fitting function takes 
into account mistag dilution and is convoluted with the 
At resolution function described above. We set St and 
Ct to zero since the transverse polarisation in the fit is 
small. 

We perform an unbinned extended maximum likeli- 
hood (ML) fit that assumes the event types mentioned 
previously. The results of the fit are 246 ± 29 signal 
events with f L = 0.99 ± 0.03, S L = -0.42 ± 0.42 and 
Cl = —0.17 ± 0.27. There is a bias on the yield coming 
from the neglect of correlations in the fit (six events) and 
B-background modeling (16 events). The former is esti- 
mated using MC simulations and the latter is dominated 
by the change in signal yield when the a\p component is 
to allowed to vary in a fit to the data. The signal yield re- 
mains stable when allowing the yield of other background 
types to vary. The corrected signal yield is 224 ± 29 
events. Figure ^ shows distributions of toes > cosOi and 
to w ± w o for Lepton and Kaon 1 tagged events, enhanced in 
signal content by cuts on the signal-to-background likeli- 
hood ratios of the discriminating variables not projected. 
The additional cuts retain 0(15%) of the signal events in 
the analysis sample. For Toes and AE we show a projec- 
tion of the data for all tag categories; in these plots we 
retain 0(60%) of the signal events in the analysis sample. 
Figure shows the raw At distribution for B° and B 
tagged events. The time-dependent decay-rate asymme- 
try A = (R(At)-R(At))/(R(At) + R(At)) is also shown, 
where R(R) is the decay-rate for B° (B°) tagged events. 

The nominal fit does not account for non-resonant 
background. If we add a non-resonant component of 
B — ► pirir events to the likelihood, the fitted signal 
yield changes by less than 11% (90% C.L.). Any pos- 
sible B — > 47r component would be significantly smaller. 
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FIG. 1: The distributions for a sample of events enriched 
in signal for the variables (a) tties, (b) AE, (c) cosine of 
the p helicity angle and (d) m w ± w o. The dotted line is the 
projection of the sum of backgrounds and the solid line is 
the projection of the full likelihood. For tties we show the 
projections for (top line) all and (bottom line) Lepton and 
Kaon 1 tagged events. 



The dominant systematic uncertainties in the yield arise 
from the assumed B-background branching fractions (20 
events) and the uncertainty on the fraction of SCF events 
(14 events). The uncertainty on the estimated fraction 
of misreconstructed events is extrapolated from a con- 
trol sample of fully reconstructed B° — > D~ p + decays. 
A 10% systematic error on the branching fraction comes 
from 7T° reconstruction. The dominant systematic er- 
ror on Jl is from the uncertainty in PDF parameteriza- 
tion (±0.03). We vary CP- violation in the B background 
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FIG. 2: The At distribution for a sample of events enriched in 
signal for (a) B° and (b) B° tagged events. The dashed line 
represents the sum of backgrounds and the solid line repre- 
sents the sum of signal and backgrounds. The time-dependent 
CP asymmetry A (see text) is shown in (c), where the curve 
represents the asymmetry. 



within reasonable limits. This is the main systematic un- 
certainty on the CP results: 0.08 (0.11) on S L (C L ). The 
systematic uncertainty on Sl (Cl) from B-background 
branching fractions is 0.02 (0.03). Uncertainty in the 
vertex-detector alignment contributes an error of 0.06 
(0.04) on S L (C L ). In half of the SCF events the mis- 
reconstructed signal contains at least one wrong track; 
the difference in resolution function for these events cor- 
responds to an uncertainty of 0.03 (0.01) on Sl (Cl). 
The uncertainty in the parametrization of the likelihood 
contributes an error of 0.05 (0.02) on Sl (Cl)- We esti- 
mate the systematic error from ignoring interference with 
non-resonant modes and a\n to be 0.02 on Sl and Cl-, 
and 2.4% on the signal yield. The uncertainty from pos- 
sible CP-violation in the doubly-Cabibbo-suppressed de- 
cays on the tag side of the event is assumed to be 
the same as for B°(B°) -► tt+tt": 0.012 (0.037) for S L 
(Cl)- We also apply a correction to account for possible 
dilution from B-background (5%) and SCF (3%) events. 
Our results are 



BR(B° 



P + p-) 

h 

C L 
Sl 



(33 ± 4(stat) ± 5(syst)) x 10~ 6 
0.99 ± 0.03(stat) ±{j;^(syst) 
-0.17 ± 0.27(stat) ± 0.14(syst) 
-0.42 ± 0.42(stat) ± 0.14(syst) 



The correlation coefficient between Sl and Cl is —0.016. 
We average this branching fraction with the less precise 
result from Ref. LJ, taking into account correlations 
where appropriate [13j , to obtain the final value of (30 ± 
4±5) x 10 -6 . This measurement supersedes the previous 
BABAR result presented in Ref. Q . 

Using the Grossman-Quinn bound Q, 03 with the 
recent results on B — ► p ± p°,p°p from |3| we limit 
\a c g — a\ < 13° (68% C.L.). Ignoring possible non- 
resonant contributions, and I — 1 amplitudes one 
can relate the CP parameters Sl and Cl to a, up to a 
four-fold ambiguity. If we select the solution closest to 
the CKM best fit central value of a = 95 - 98° 0, the 
measured CP parameters of the longitudinal polarization 
correspond to a = 102l 1 2(stat)^(syst) ± 13(penguin)°. 
Figure 01 shows the confidence level as a function of 
ateff = arcsin(S , L/^/l — C|)/2 for this result, (dotted) 
taking into account systematic uncertainties and (solid) 
also including the penguin contribution. We exclude val- 
ues of a between 19° and 71° (90% C.L.). 
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confidence level for this result. 
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